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Abstract The origin of hydrodynamic turbulence in rotating shear flow is a long standing 
puzzle. Resolving it is especially important in astrophysics when the flow angular momen- 
tum profile is Keplerian which forms an accretion disk having negligible molecular viscosity. 
Hence, any viscosity in such systems must be due to turbulence, arguably governed by mag- 
netorotational instability especially when temperature T> 10^. However, such disks around 
quiescent cataclysmic variables, protoplanetary and star-forming disks, the outer regions of 
disks in active galactic nuclei are practically neutral in charge because of their low temper- 
ature, and thus expected not to be coupled with the magnetic field appropriately to generate 
any transport due to the magnetorotational instability. This flow is similar to plane Couette 
flow including the Coriolis force, at least locally. What drives their turbulence and then trans- 
port, when such flows do not exhibit any unstable mode under linear hydrodynamic pertur- 
bation? We demonstrate that the threedimensional secondary disturbance to the primarily 
perturbed flow triggering elliptical instability may generate significant turbulent viscosity 
ranging 0.0001 < i/(< 0.1 to explain transport in accretion flows. 
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1 INTRODUCTION 

One of the main problems behind the origin of hydrodynamic turbulence in shear flow is that there is a sig- 
nificant mismatch between the predictions of linear theory and experimental data. For example, in the case 
of plane Couette flow, laboratory experiments and numerical simulations show that the flow may be turbu- 
lent at a Reynolds number as low as Re ^ 350, while according to the linear theory the flow should be stable 
for all Re. Similar mismatch between theoretical results and observations is found in astrophysical con- 
texts, where the accretion flow of neutral gas with Keplerian angular momentum profile, which essentially 
behaves like rotating shear flow, is a common subject. Examples of such flow systems are accretion disks 
around quiescent cataclysmic variables dGammie & Menou 19981 1, protoplanetary and star-forming disks 
(IBlaes & Balbus 1994| l, and the outer regions of disks in active galactic nuclei (Me nou & Quataert 200 l| l. 

A Keplerian accretion disk flow having a very low molecular viscosity must generate turbulence and 
successively diffusive viscosity, which support the transfer of mass inwards and angular momentum out- 
wards. However, theoretically this flow, in absence of magnetic field, never exhibits any unstable mode 
which could trigger turbulence in the system. On the other hand, the laboratory experiments of Taylor- 
Couette systems, which are similar to Keplerian disks, seem to indicate that although the Coriolis force 
delays the onset of turbulence, the flow is ultimately unstable to turbulence for Reynolds numbers larger 
than a few thousand (Richard & Zahn 2 001) , even for subcritical systems. Indeed, Bech & Anderson (1997) 
see turbulence persisting in numerical simulations of subcritical rotating flows for large enough Reynolds 
numbers. 

How does shearing flow that is linearly stable to perturbations switch to a turbulent state? Since last 
decade, many authors including ourselves have come forward with a possible explanation of this fact 
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based on bypass transition (see, IButler & Farrell 1992', 'Reddy & Henningso n 1993] ITrefethen et al. 19^31 
Chagelishvili et al. 2003 Umurhan & Regev 2004 Mukhopadhyay et al. 2005J and references therein) 



where the decaying linear modes show an arbitrarily large transient energy growth at a suitably tuned 
perturbation. In lieu of linear instabilities e.g. magnetorotational instability, the transient energy growth, 
supplemented by a non-linear feedback process to repopulate the growing disturbance, could plausibly sus- 
tain turbulence for large enough Reynolds numbers. 

The behavior of shear flows, however, in the presence of rotation is enormously different compared 
to that in absence of rotation. The Coriolis effect is the main culprit behind this change in behavior 
killing any growth of energy even of transient kind in the presence of rotation. In the case of shear flow 
with a varying angular velocity profile, e.g. Keplerian accretion flow, the above mentioned transient en- 
ergy growth is insignificant for threedimensional perturbations. To overcome this limitation, it is nec- 
essary to invoke additional effects. Various kinds of secondary instability, such as the elliptical insta- 
bility, are widely discussed as a possible route to self-sustained turbulence in Unearly perturbed shear 
flows (see, e.g. IRerrehumbert 1986 Bayly 1986 Craik & Criminale 1986, Landman & Saff man 1987) 



Hellberg & Orszag 1988, Waleffe 1989. ,Craik 1989. ,Le Dizes et al. 1996. .Kerswell 2002,) . These effects, 



which generate threedimensional instabilities of a twodimensional flow with elUptical streamlines, have 
been proposed as generic mechanism for the breakdown of many twodimensional high Reynolds num- 
ber flows whose vortex structures can be locally seen as elliptical streamlines. Recently, one of the 
present authors has studied the secondary perturbation and corresponding elliptical vortex effects in ac- 
cretion disks and pinpointed that they can be the seed of threedimensional hydrodynamic instability 
( [Mukhopadhyay 2006 ) . Subsequently, by numerical simulation, this has been shown to be one of the pos- 
sible sources to generate turbulence to form large objects from the dusty gas surrounding a young star 
jCuzzi 20071 lOrmel et al. 2008 1). Moreover, vortex generation in the unmagnetized protoplanetary disks has 
been furnished by hydrodynamic turbulence dde Val-Borro et al. 2007 ) which leads to planet formation, and 
angular momentum transport in disks. However, whether they lead to non-linear feedback and threedimen- 
sional turbulence are yet to be shown explicitly. 

Here we plan to show in detail that threedimensional secondary perturbation generating large growth 
in the flow time scale may generate significant turbulent viscosity in rotating shear flows, more precisely in 
plane shear flows with the Coriolis force. The plane shear flow with the Coriolis force essentially behaves 
as a local patch of a rotating shear flow. Possibility of significant turbulent transport in such flows by three- 
dimensional perturbation opens a new window to explain accretion process in flows which are neutral in 
charge. In particular, we address the issue of deriving turbulent viscosity and the Shakura-Sunyaev viscosity 
parameter a ( Shakura & Syunyaev 1973 1 from a pure hydrodynamical perspective Q. This is important for 
understanding accretion flows in cold charge neutral medium. 

It is important to note that transition to turbulence is not a unique process, but it depends on the initial 
condition/disturbance and the nature of the flow ( jSchmid & Henningson 20()T ICriminale et al. 20031 ). In 
fact, it is known that even in the presence of secondary instability, linearly unstable base flows may reach 
to a non-turbulent saturated state. However, turbulence definitely belongs to the nonlinear regime and it is 
exhibited only in the situations when large growth of perturbation switches the system over the non-linear 
regime. As our present goal is to understand the possible origin of hydrodynamic turbulence, we consider 
those situations when large energy growth governs non-linearity. 

The paper is organized as follows. In the next section, we first recall the perturbation established pre- 
viously (Mukhopadhyay 2006) due to secondary disturbance in the Keplerian flow and then discuss the 
range of corresponding Reynolds number and the solutions. Subsequently, we estimate the corresponding 
turbulent viscosity of hydrodynamic origin in §3. We end in §4 by discussing implications of our results. 



2 PERTURBATION AND RANGE OF REYNOLDS NUMBER 

Considering a twodimensional velocity perturbation w = {wx{x, y, z, t), Wy{x, y, z, t), 0), and pressure 
perturbation Pp{x, y, z, t) in a small section of the Keplerian shear flow/disk, the linearized Navier-Stokes 
and continuity equations for the incompressible fluid with plane background shear in the presence of a 

' A preliminary calculation of such a has been appeared in a collected volume of Gravity Research Foundation 
^Mukhopadhyay 2008) . 
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Coriolis component can be written in dimensionless units as (see |Mukhopadhyay et al. 2005| for a detailed 
description) 

— = 2nw - ^ + — 

dt ^ dx Re 



= 2nwy - ^ + — V^u-^, (1) 



dwx dwy 

1 

dx dy 



= 0. (3) 



We consider the standard no-slip boundary condition such that Wx = Wy = at a; = ±1 and according 
to the choice of variables in the coordinate system 57 = 1/q. Here {x, y, z) is a local Cartesian coordinate 



system centered at a point (r, (f) in the disk ( Mukhopadhyay et al. 2005 1 such that dr ~ x and rd(f> — y 



When the Reynolds number is very large, the solution of eqns. ([T]i, © and ^ are given by 
( [Mukhopadhyay et al. 2005| l 



k k 
Wx = C-^ s\a{kxX + kyy), Wy = sin(fc^x + kyy) (4) 



where ( is the amplitude of vorticity perturbation, kx and ky are the components of primary perturbation 



wavevector and I — ^Jk^ + ky. Under this primary perturbation, the flow velocity and pressure modify to 

U = U'' + w = {wx,~x + Wy,0) ^ A.d, P^p + pp, (5) 

where U^, p are background velocity and pressure respectively, A is a tensor of rank 2. Here kx — kxo+kyt, 
which basically is the radial component of primary perturbation wavevector, varying from — oo to a small 
number, where kxo is a large negative number: |fca;o| ^ i?e^/^ ^ tmax ( [Mukhopadhyay et al. 2005| l. 

Now we concentrate on a further small patch of the primarily perturbed flow such that the spatial scale 
is very small compared to the wavelength of primary perturbation satisfying sin(A;j;2; + fcj,?/) ~ kxX = f< 1. 
In fact, / 1 at close to the boundary of the patch when i/ — > and 2'K/ky, and at an intermediate location 
/ ^ 1. As |fca; I varies from a large number to close to unity, the size of the primary perturbation box in the 
cc-direction is l/kx< 1 when ky ~ 1, fixed. Hence, this further small patch must be confined to a region: 
— a< x< a, when f /\kxo\< a< f. Clearly, in this patch, U in eqn. (O describes a flow having generalized 
elliptical streamlines with e = {kx/l)"^, a parameter related to the measure of eccentricity 0, running from 
to 1 as the perturbation evolves. It was already shown ( [Mukhopadhyay 2006| l that a secondary perturbation 
in this background may grow exponentially leading the flow unstable. We use this unstable flow in §3, which 
was extensively discussed earlier ( Mukhopadhyay 2006 ), to derive i^t and a. 

As we focus on the secondary perturbation at a small patch of the primarily perturbed shearing box, the 
variation of primary perturbation appears insignificant in the patch compared to that of the secondary one. 
Depending on the primary perturbation wavevector at a particular instant, the size of the secondary patch 
is appropriately adjusted. In fact e varies very very slowly and marginally deviates from unity in the time 
interval when kx varies from kxo (large negative) to, say, —10. Even when kx tends to —3, e changes to 
~ 0.9 only. Therefore, e and thus A practically remains constant. 

2.1 Range of Reynolds number 

Due to consecutive choice of small boxes/patches, the Reynolds number in the secondary flow is restricted 
with a particular choice of that in the primary flow. Here in the interest of clarity, we work with the original 
dimensioned units. The Reynolds number at the primary box is defined as 

UoL qnpL^ 

Rep = = , (6) 



^ Note that e is a parameter related to the measure of eccentricity but not the eccentricity itself. 



4 



Mukhopadhyay & Saha 



where 2L is the box size in the x-direction and 2Uo is the relative velocity of the fluid elements in the 
box between two walls along the y-direction. Now we recall the secondary perturbation at a smaller patch, 
extended from —Lg to +Ls, such that |is | aL. To meet our requirement fim{kxX + ky-y) ~ k^x + ky-y, 
we remind that the small patch size needs to be adjusted. Therefore, the Reynolds number at the secondary 
box is given by 

R,^^s^^ i^h^, (7) 

Hence, 

Re, ^ a^^ p- 

At the beginning of the primary perturbation k^ = fcxo and thus e = 1. At this stage, the secondary box 
size Lg — Lf /kxo and Rep> k^^Reg- With time k^ decreases in magnitude but e deviates little from 
unity until k^ ^ —3 when e — 0.9. Hence A can be considered constant approximately as described 
above. At this stage Rcp > 9 Res, atleast an order of magnitude higher than Re,- If the energy growth 
due to primary perturbation is maximized for k^ = kx,min = tt (Muk hopadhyay et al. 2005) l, then the 
range of Re for the secondary perturbation is given by Rep /^/fc^o~ Rgs< Rcp p /lO. At kx — tt, Rcs 
is atleast an order of magnitude lower than Rep. When kx,mim — 1, Rep ^ Re, for / ~ 1. In general 
-Sep P/kxo< Res< Rep P /k"^ ,^^^. 

2.2 Solution 



Following previous work ( [Mukhopadhyay 2006 1, the general solution for the evolution of secondary per- 



turbation in the flow discussed above can be written in terms of Floquet modes 

Ui{t) ^ ex.p{at) ft{4})cxp[i{kix + k2y + ksz)], (9) 

where (j> ~ wt, fi{4>) is a periodic function having time-period T — 2tt/w, a is the Floquet exponent, 
ki,k2, k^ are the components of wavevector of the secondary perturbation. Note that a is different at differ- 
ent e. Clearly, if cr is positive, then the system is unstable. The detailed solutions were discussed elsewhere 



( Mukhopadhyay 2006 ) what we will not repeat here. 



In principle, kx varies with time and thus A does so. Thus, generalizing the solution (|9]) for a (slowly) 
varying A, we obtain 

Ui{t) = exp (^J a{t) dt^ fi{(j)) exp[i(kix + k2y + k-^z)], (10) 

where (f> — J ^{t) dt. The eqns. (|9]l and ( fTOl i practically describe the solutions for the entire parameter 
regime exhibiting elliptical vortices which are very favorable for the elliptical instability to trigger 

For the present purpose, the physically interesting quantity is the energy growth of perturbation which 
is given by 



where Y,{t) — J a{t) dt and t = {kx — kxa)/ky. As kx{t) varies from a large negative value, kxo, to 0, 
t increases from to t^ax — —kxo/ky. Thus, the energy growth is controlled by the quantity as 
fi{(l>)/fi{0) simply appears to be a phase factor Therefore, our aim should be to evaluate E for various 
possible perturbations. 

Let us specifically concentrate on the Keplerian accretion flows. Figure[T^ shows the variation of max- 
imum velocity growth rate, (Tmax, as a function of eccentricity parameter, e, for the various choices of 
amplitude of vorticity, (. By "maximum" we refer the quantity obtained by maximizing over the vertical 
component of the wavevector, k^. At large e (as well as large kx), when C is large, the background flow 
structure. A, is elliptical with high eccentricity. Therefore a vertical perturbation triggers the best growing 
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Fig. 1 (a) Variation of maximum velocity growth rate as a function of eccentricity parameter. 
Solid, dotted, dashed and long-dashed curves indicate the results for C — 0.01,0.05,0.1,0.2 
respectively (i Mukhopadhyay 2006| l. (b) Variation of E as a function of time for fc^-o = —10^, 
when various curves are same as of (a), (c) Same as (b) but for kxo — —10"*. Other parameters 

are kyo — 1, fcio = 0, |feo| = 1, and q ~ 3/2. 
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mode into the system. However, with the decrease of (, A approaches to that of the plane shear and thus the 
growth rate decreases significantly. At this stage, the corresponding best perturbation is threedimensional 
but not the vertical one. 

At small e (and then small fc^), when ( is large the eccentricity of the background elliptical flow de- 
creases significantly, and thus the growth rate decreases. In this low eccentric flow, the best growth rate 
arises due to the twodimensional perturbation. On the other hand, when is small, the background reduces 
to that of the plane shear flow. Therefore, the growth rate increases according to the shearing effects, as 
described by [Mukhopadhyay et al. 2005 An interesting fact to note is that except the case of small e (kx) 



with a large the growth rate maximizes for the threedimensional perturbation. Moreover, at a large ^ and 
a large e, the best growth rate arises due to a vertical (or almost vertical) perturbation. 

As the accretion time scale is an important factor, for the present purpose, physically interesting quantity 
is S rather than a itself. Figures [T]3,c show the variation of S as a function of t at various C,. As the 
perturbation evolves with time, the corresponding E increases. It is also clear that S and then corresponding 
growth increases with the increase of \kxo \ (and then Re), i.e. the increase of accretion time scale, in addition 
to the increase of (. In Table 1, we enlist the approximate values of maximum growth factor, as follows 
from eqn. ( fTTT i. corresponding to Emax ~ /o for 'he cases shown in Figs.[T]3,c. When k^o = —10^, 

Rcp ~ 10^^ (as Rep ^ tf-^^ax ^ k^a) ^i^^ from eqn. dHJ Res{f = 1)> 10^, the maximum growth factor is 
significant for a large amplitude of vorticity perturbation i.e. C, > 0.1. However, the growth factor increases 
with the increase of Rep and when Rep ~ 10^^, and then Res{f = 1)> 10^, it is quite significant for an 
amplitude of vorticity perturbations as small as 0.05. Therefore, it appears that a suitable threedimensional 
secondary perturbation efficiently triggers elliptical instability and possible turbulence in rotating shear 
flows including accretion disks. 

Table 1 

Maximum energy growth corresponding to cases shown in Figs.[T]D,c 



\kxo\ 


c 


^max 


G^max 


10^ 


0.2 


6.1 


2 X lO'' 


10^ 


0.1 


5.2 


3.3 X 10^ 


10^ 


0.05 


4.43 


7 X 10^ 


10^ 


0.01 


1.97 


52 


10* 


0.2 


3.65 


1500 


10-* 


0.1 


3 


400 


10^ 


0.05 


2.9 


330 


10* 


0.01 


1.27 


13 



3 TURBULENT VISCOSITY 

Here we attempt to quantify the turbulence by parametrizing it in terms of the viscosity. This is essentially 
important, as explained in §1, in flows like astrophysical accretion disks, where molecular viscosity is 
negligible, to explain any transport therein. 

The tangential stress at a point (r, 0) of a rotating flow exhibiting turbulence is 

dn 

WrA = vt r— = -vt q^, (12) 
dr 

where vt is the turbulent viscosity and ft — fioC'V'^o) Note that 5 = 3/2 for the Keplerian angular 
velocity profile. The perturbation described above is expected to govern the nonlinearity after certain time, 
say tg. We also assume that the nonlinearity leads to turbulence attributing the fact that at the initiation 
of turbulence the eddy velocity is same as the perturbation velocity. Therefore, we obtain the averaged 
tangential stress due to perturbation alt ^tg 

Ux{tg)uy{tg)dxdy, (13) 



4ttL, 



-is 
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where we remind that the azimuthal flow is considered to be periodic iny — 27r/fc2. 
Now combining eqns. (fT2b . ( fTSl ) and after some algebra we obtain 

= (14) 

where Txy = J Wxy dxdy, M = Vtx/cs and vt denotes the averaged vt in the small section, computed here 

?At — tg. 

Without any proper knowledge of turbulence in Keplerian flows which arise in accretion disks, Shakura 
& Sunyaev ( [Shakura & Syunyaev 1973| parametrized it by a constant a considering Wrtp to be proportional 
to the sound speed, c^, given by 

Wr4, = -acl (15) 

a is called the Shakura-Sunyaev viscosity parameter They assumed that the small section under considera- 
tion to be isotropic so that scaled the characteristic length It of turbulence in terms of the largest macroscopic 
length scale of the disk, i.e. half-thickness h, and the eddy velocity of turbulence vt in terms of sound speed 
Cs- Thus they defined the turbulent viscosity 

Vt = = acs/i, (16) 

where It — aih, vt — ayCs, a — aiav/3- Obviously ai < 1. If the turbulent velocity becomes supersonic, 
then shock forms and reduces the velocity below the sound velocity which assures a^, < 1. Therefore, a< 1. 
From eqns. ( fT4] i and ( fTSI l we write 

T-y (17) 



where a denotes the averaged a in the small section. Therefore, if we know the structure of the flow, then 
we can compute the turbulent viscosity due to various perturbations. As we consider the size of the section 
to be very small, a and Dt are effectively equivalent to a and i^t at a particular position in the disk. Below we 
compute Txy for the various secondary perturbations and the corresponding turbulent viscosities, at least in 
certain approximations. 

3.1 Secondary perturbation evolves much rapidly than the primary one 

From eqn. (|9]l we can write the velocity perturbation components 

Ux{x,y) ^j;e'^*/x(0)sin(fcia; + fc2j/ + /c3 2), 

Uy{x,y) = Aye''*'fy{(j))siii{kix + k2y + k3z), (18) 

where Ax and Ay are the amplitudes of perturbation modes, A;io,^20 are the radial and the azimuthal 
components respectively of the secondary perturbation wavevector at i = 0, Ax and Ay can be evaluated 
by the condition that the velocity components of the secondary perturbation reduce to that of the primary 
perturbation att — (at the beginning of the evolution of secondary perturbation) given by 

_ s\n{kx{e)x + kyy) 
sm(fcioa; + fc2oy + fc30^) 



where kxic) = — e)ky, C is of the order of unity (for details see Mukhopadhyay et al. 2005 

Mukhopadhyay 2006[ ). Therefore, from eqn. ( fT3] l 
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Now by considering a typical case with ky = 0.71, vt and a can be computed as functions e (kx), when we 
know the time of evolution of the secondary perturbation tg. 

Figure |2]describes vt and a according to eqns. (fl4] l. ( [TtI i and (l20l i for various disk parameters. As the 
primary perturbation evolves, elliptical vortices form into the shearing flow which generate the turbulent 
viscosity under a further perturbation. Figure shows that the viscosity varies with the eccentricity of 
vortices. At a very early stage when the primary perturbation is effectively a radial wave and e — > 1, the 
maximum velocity growth rate due to secondary perturbation, a„iax (shown in Fig.[T^), and the correspond- 
ing turbulent viscosity are very small, independent of the value of C,. With time, the primary perturbation 
wavefronts are straightened out by the shear until t ~ tmax, when the perturbation becomes effectively an 
azimuthal wave and e — s- 0. At this stage, amax and the turbulent viscosity due to the secondary perturba- 
tion become zero again. This feature is clearly understood from eqn. (l20l i. However, at an intermediate time 
when kx{e) is finite, vt may be ^ 0.005 even in a moderately slim disk with h{r)/r ~ 0.05, when the time 
of evolution of secondary perturbation tg = 10. This tg is considered to be the time at which turbulence 
is triggered in the system. Figures |2j5-d show the variation of h't and a with the eccentricity of vortices at 
various C, when tg = 10, 100. It is interesting to note, particularly for tg = 100, that with the increase of C, 
first viscosity increases then decreases. This is understood from the underlying energy growth rate shown 
in Fig.[T^, when the readers are reminded that a = (t(C, e). Note that the qualitative behavior of i^t is same 
as that of a. If we look at a typical case with = 0.05 where a = amax at e = 0.86 which corresponds to 
kx = —1.76, then a and I't computed att — tg are for Res< Rcp ^ 10^. 

3.2 Secondary perturbation over the slowly varying primary perturbation 

In principle, the primary perturbation may vary with time during the evolution of secondary perturbation. 
By numerical solutions, simultaneous evolution of the primary and the secondary perturbation along with 
the corresponding energy growth has already been discussed earlier (Mukhopadhyay 2006). For the conve- 
nience of analytical computation of viscosity, here we consider the regime of slow variation of the primary 
perturbation compared to the secondary one. Hence we recall eqn. ( fTOl i and write the velocity perturbation 
components 

Ux Uxj:{x,y) = Bxe^'-*^ fxi(l))sin{kix + k2y + ks z), 

Uy Uy^{x,y)=Bye^^*^fy{(l))sm{kix + k2y + k3z), (21) 

with (j) = J m{t)dt. The amplitudes of perturbation modes Bx and By can be evaluated by the initial 
condition of secondary perturbation. The secondary perturbation could trigger elliptical instability only after 
significant vortex forms in the flow due to the evolution of primary one. At the beginning of the evolution 
of primary perturbation kxo — oo (we choose the cases kxo = —10^ and — 10"*) which corresponds 
to e — > 1 and thus effectively a plane shear background when ( is small (see jMukhopadhyay 2006 1. In 



absence of vortex, this can not trigger elliptical instability under a secondary perturbation. As kxo decreases 
in magnitude, e deviates from unity giving rise to a background consisting of elliptical vortices. Above 
certain e = ec, the secondary perturbation does not have any effect to the primarily perturbed flow and 
Uxj^ and Uy^ reduce to the primary perturbation. We hypothesize that ec ~ 0.9999. Hence, Bx and By are 
computed in a similar fashion as in §3. A given by 

Bx = Cirrh^TT^T, By = -C- 



(J ^ &vcL{kx{ec) X + ky y) ^22^ 
sin(fcio x + k2Qy + k^o z) ' 



Hence, from eqn. ( flj] ) the stress tensor 



T (f ^ ^ .2 ^x{e-c)ky 2S„ax p 

^xyy'"caai^) ^ 2/'*(e ) 

n _ r^2 fx{4>)fy{4>) 
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Fig. 2 This is for the perturbation described in §3. A. (a) Variation of vt (dotted curve) and 
a (soHd curve) as functions of e for ( = 0.05 case described in Fig. [T^, when h{r)/r = 
0.01,0.05,0.1 respectively for the top, middle, bottom curves of a; r = 30, ky = 0.71, 
tg = 10. (b) Variation of i^t as a function of e for the cases described in Fig. [T^ with 
h{r)/r = 0.05, tg = 10, ky = 0.71, when solid, dotted, dashed, long-dashed curves correspond 
to C = 0.01, 0.05, 0.1, 0.2 respectively with |fco| = 1- (c) Same as in (b) except a is plotted in 
place of vt- (d) Same as in (c) except tg = 100. 
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Fig. 3 This is for the perturbation described in §3.B. Variation of vt (dotted curve) and a (solid 
curve) as functions of h{r) j r for cases shown in Figs. ^p,c, when the curves from top to bottom 
correspond to C = 0.2, 0.1, 0.05, 0.01 with r = 30 for (a) = -10^ (b) = -10''. Other 
parameters are ky = 1, — 0.9999. 



where kx reduces to zero sit — tmax, which corresponds to the beginning of turbulence when S = Smax- 
It is found from Fig. |3]that in a thin disk with h{r)/r = 0.01, a at r = 30 may be as high as > 0.1 
for kxQ = — lO'^ when C, is very large. Although the viscosity decreases with the decrease of (, a still 
may be ~ 0.001 when ( ~ 0.05. The turbulent viscosity decreases in a considerably thicker disk, but still 
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a ~ 0.003 at h{r)/r = 0.1 when Q = 0.2. For C > 0.1, vt> 0.001 when = -10^ The values of vt 
and a both decrease when jfcxol decreases to 10^, which is expected from Table 1 as well. In this case, a 
significant turbulent viscosity generates only at a large C, = 0.2. 

4 IMPLICATIONS AND DISCUSSIONS 

Above results verify that at a range of e, the threedimensional growth rate due to secondary perturbation 
in rotating shear flow of the Keplerian kind is always real and positive and corresponding growth may be 
exponential and significant enough, at least for a suitable choice of C, and/or Re, to trigger non-linearity 
and then plausible turbulence in the flow time scale. With the increase of k^o ('^ Rcp^^'^), the effect due to 
elliptical instability increases, and thus corresponding growth does so. 

As this growth is the result of threedimensional perturbation, underlying perturbation effect should sur- 
vive even in the presence of viscosity. There are many important natural phenomena where the Reynolds 
number is very large. In astrophysical accretion disks, what applications are essentially considered in the 
present paper. Re always could be > lO^'' because of their very low molecular viscosity. Therefore, the 
present mechanism is certainly applicable to such disk flows to resolve their turbulence puzzle when es- 
pecially it is cold and neutral in charge and thus not a very plausible candidate for the magnetorotational 
instability. On the other hand, we suggest that the subcritical transition to turbulence in Couette flow may 
be the result of secondary perturbation which triggers elliptical instability modes into the system. 

We have tried to estimate the corresponding hydrodynamic turbulent viscosity. We have aimed to quan- 
tify the amount of turbulence through this using the perturbations as the source of turbulence. We report 
here an observable range of viscosity obtained for the typical thin accretion disks and with reasonable val- 
ues of flow vorticity. In place of r = 30, if we choose the shearing box at a large distance from the central 
object, say at r = 500, then the computed a naturally decreases three orders of magnitude [see eqn. (fT7]l1. 
We show by an extensive analysis the dependence of viscosity on the aspect ratio [h/r) of the flow. The 
values of vt and a increase quite rapidly as the disk becomes thin to thinner From eqns. (fT4l i and (fTTb and 
with the results given in Figs. |2] and [3] we find that it still might be as large as 10^^ for a thin disk even at a 
large distance, say, r — 500. 

While some earlier laboratory experiments (e.g.l Richard & Zahn 20011 ) predicted sub-critical transition 
to turbulence and then transport in hydrodynamical shear flows like accretion disks, experiments by Ji et 
al. (2006) have argued against it. Non-detection of turbulence and then any angular momentum transport 
of purely hydrodynamic origin could be due to the following facts. Maximum Reynolds number in this 
experiment is 2 x 10^ whereas the cold disks such as the protoplanetary disks have Reynolds number 
^ 10^^. However, the critical Reynolds number for these systems could be ^ 10^ — 10^ or more. It can 
be easily understood with a very simple example that as Re increases, the amplitude of vortices increases 
which are indeed clear from the Figs. 7 and 8 given by Mukhopadhyay et al. (2005). Let us consider a 
2D perturbation in an inviscid incompressible flow where the vorticity V x v is exactly conserved, when 
V — ivx + jvy. Therefore, at f = t,nax — tg, when the perturbation growth is maximum at t = t^ax, 
the amplitude of vorticity ^ ^ ^ Re^/^. As ut and a are directly proportional to C^, they scale as 
Re^l^ at t = tmax — tg- Therefore, if Re decreases three orders of magnitude, then vt decreases in two 
orders. Moreover, the perturbation stabilizes at a thicker disk. Indeed we find that the viscosity decreases, 
as h{r)/r increases. Dimension of confined liquid in the experiments by Ji et al. (2006) may not be typical 
of astrophysical disks or rings, when they may have a large aspect ratio ^ 2, whereas the astrophysical 
disks and ring systems are normally thin (with aspect ratio < 1). Obviously a huge gap exists between 
experiments and the real observations. 

By numerical simulations, the formation and evolution of vortices in a hydrodynamic shearing-sheet 
have already been studied by Johnson & Gammie (2005) and they suggested it to be a possible mechanism 
for angular momentum transport in low-ionization disks at high resolution. It has been argued that there 
must be a mechanism to inject vorticities into the disk, and the vortices must not decay rapidly due to three- 
dimensional instabilities, to sustain the transport. We show that the vortices may sustain in threedimension 
at least in the time scale of interest, where this is applicable for accretion disks. Indeed, Cuzzi and his collab- 
orators jCuzzi 2 007 , Orm el et al. 20 08 ) have argued, by numerical simulations, that the elliptical instability 
may lead to turbulence to from the dusty gas surrounding a young star. Also the vortex generation and 
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then the angular momentum transport has been shown to occur in the unmagnetized protoplanetary disks 
dde Val-Borro et aL2 007 ) by hydrodynamic turbulence. However, other simulations jShen et £17 2006 ) do 
not find significant transport. The nonoccurence of significant transport in simulations, in our view, is due 
to lack of resolution needed to capture the turbulence. Indeed, the later authors have mentioned that for their 
calculations it is difficult to define an effective Reynolds number, since the numerical dissipation is a steep 
function of resolution. With a particular non-linear solution, Balbus & Hawley (Balbus & Hawley 20061 
have shown that perturbation decays asymptotically. They also have argued that as the nonUnear term in 
the equation for the incompressible flow itself vanishes expUcitly, the solution can not lead to nonlinearity 
and then turbulence. However, this does not guaranty that every solution does so. They themselves have 
also mentioned that secondary instabilities may still spoil their conclusion. Indeed the coupling between 
the secondary and primary modes was shown earlier not to allow the nonUnear term to vanish resulting in a 
possible nonlinear transition to turbulence ( |Mukhopadhyay 2006 1. 

It is interesting to note that the modal instability via the bypass mechanism (and then with a secondary 
perturbation superimposed) arises in these systems from a subtle interplay of the non-normality of the 
perturbation modes and the non-linearity of the Navier-Stokes equation and this in turn gives rise to the 
turbulence in the system. As the turbulence and corresponding transport is inevitable in these systems, the 
corresponding a may not be just inversely proportional to the critical Reynolds number (as predicted earlier 
( Lesur & Longaretti 2005) ). Previous theoretical studies (Mukhopadhyay et al. 200 5 ) have shown that the 
Keplerian flow may render a transition to the turbulent regime at a Reynolds number ^ 10^ and turbu- 
lence might have just started at this critical Reynolds number It is to be seen now whether all shear flows, 
exhibiting subcritical turbulence in the laboratory, do exhibit large growth due to secondary perturbation. 
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